Periodic revival of entanglement of two strongly driven qubits in a dissipative cavity 
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We study the dynamics and decoherence of a system of two strongly driven qubits in a dissipative 
cavity. The two qubits have no direct interaction and are individually off-resonantly coupled to 
a single mode of quantized radiation. We derive analytical solutions to the Lindblad-type master 
equation and study the evolution of the entanglement of this system. We show that with non-zero 
detuning between the quantum and classical fields, the initial decay of the entanglement is followed 
by its revival periodic in time. We show that different Bell states follow evolutions with different 
rates. 
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Entanglement is a resource essential for successful im- 
plementation of quantum information processing [1], and 
engineering sustainable entanglement is a necessary re- 
quirement for any physical realization of a quantum com- 
puter. Since a number of solid-state qubits have been suc- 
cessfully realized and single qubit operations have been 
demonstrated [2], the attention of experimentalists is now 
mainly focused on interaction between qubits. It is there- 
fore of vital importantce to study processes which can 
generate or alter entanglement through an interaction 
between two or more qubits. 

In many realizations it is easier to couple solid-state 
qubits via an optical or a microwave cavity rather than 
directly. Entanglement between qubits coupled to optical 
cavities has been studied in a very broad context. It is 
known that dissipation, both in qubits and in the cavity, 
can very quickly lead to disentanglement of two qubits 
[3-5]. This process, known as the entanglement sudden 
death (ESD), was confirmed experimentally [6, 7]. For 
carefully chosen initial conditions and system parame- 
ters, ESD can be followed by the entanglement revival - 
entanglement sudden birth [8-10], which, however, was 
found to never reach the initial (before ESD) degree of 
entanglement. Later it was shown that two directly in- 
teracting qubits can become entangled via a spontaneous 
decay [11]. Numerical results supporting the claim that 
under spontaneous atomic decay some atom-field entan- 
glement can be generated have been presented in Ref. 
[12]. 

In this Letter, we consider two qubits which are 
strongly coupled via a dissipative cavity, so that the di- 
rect interaction between the qubits is negligible. The 
qubits are driven by a classical ac off-resonant field. As 
a result of detuning between the classical pumping fre- 
quency and the quantum cavity eigenmode, we find a 
novel entanglement behaviour where two initially entan- 
gled qubits experience periodic entanglement drops and 
revivals. Their concurrence, despite its initial sudden de- 
crease with time, asymptotically tends to its initial value. 
Moreover, we find that once the qubits are unevenly cou- 
pled to the cavity, the decoherence free subspaces pre- 




FIG. 1: Schematic representation of the setup. 

viously discussed [13] disappear and the previously un- 
affected states also decay however now at slower decay 
rates. Our results can be applied to superconducting 
qubits coupled to a microwave cavity or to an NV center 
in diamond strongly coupled to an optical cavity. 

Below, we first present the derivation of the effective 
multi-qubit Hamiltonian in the interaction picture. Then 
we find the solutions to a single qubit interacting with a 
coherent mode of radiation in a dissipative cavity, which 
is followed by an extension of this problem to the two- 
qubit case. Finally we analyze the temporal dynamics of 
entanglement of this system. 

The model. The interaction of a qubit and a cavity is 
commonly described in terms of the Jaynes-Cummings 
Model (JCM) [14] which is one of the few interacting 
quantum systems admitting closed form solutions. JCM 
and its several variants have become a textbook tool to 
discuss coupled qubit and photon systems. Recently, it 
has been realised, that the qubit-ficld interaction with an 
additional strong driving also can be solved analytically 
[15-17] even if the cavity dissipation is also included in 
the system. Moreover, in Refs. [18-21] it has been proven 
that the solutions to the equations of motion for strongly 
driven qubits interacting through a cavity vacuum field 
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can be extended to an unlimited number of qubits, which 
can not be achieved in the simple JCM Hamiltonian. 

To keep the treatment general, we consider a Hamil- 
tonian of a system of N identical qubits coupled to a 
single-mode cavity and additionally driven by a classical 
electromagnetic field [17], 
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where Q is the level spacing of the qubits, w is the 
frequency of the eigenmode of the cavity, A and uj c 
are the amplitude and the frequency of the classical 
field, and gj is the coupling strength between the j th 
qubit and the cavity mode. In addition to that a z and 
cr ± = | (<j x ± icr 2 ') are the (linear combinations of) Pauli 
matrices and a (a*) is the annihilation (creation) oper- 
ator of the quantum field mode. Throughout the Letter 
we set h = 1. 

We assume that the qubits are driven strongly and that 
they are very stable and moderately coupled to the cavity 
mode, A » co, 5 » g » 7, where 7 stands for the qubit 
decay rates. Therefore, we can ignore the qubit dephas- 
ing or decoherence rates as well as the energy violating 
("counterrotating") Rabi Hamiltonian terms [17], a + a) 
and a~a. Additionally, we take the classical field to be 
sufficiently off- resonant, lj c 7^ w, so that wc can ignore 
the classical field-cavity coupling. 

We begin by applying an entanglement preserving 
time-local unitary transformation 

H -> W = U ] HU -iU^dtU , 

IV) -+ W) = v\<i>), 

with U = exp ^— iuj c ta)a — iuj c t cr |/2^ ■ The resulting 
Hamiltonian now takes the form 

H = H Q + Hj , 
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with A = O — lq c and 5 = uj — uj c . The interaction picture 

Hamiltonian V = e~ tH ° t Hie 1 ' H ° t upon setting the qubits 
in resonance with the classical field A = yields 
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Master equation, one qubit. Let us first focus on an 
interaction between a single qubit N = 1 and a coherent 
state of the cavity a. The evolution of this system in a 
dissipative cavity is driven by the Lindblad-type master 
equation, 
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where V (p) = 2apa) — a) dp — pa) a = 2A4 (p) —IZ(p) — 
£ (p) is the so-called dissipation operator and k repre- 
sents the cavity decay rate. 

Using the interaction picture Hamiltonian and ex- 
pressing the qubit density matrix in the |± ) basis as 
p(t) = Pkl(t)\k)(l\, k,l = ±, one can write the equa- 
tions of motion (3) for individual density matrix entries, 

p ++ (t)\a)(a\ = -ig[ae- M + ah iSt , \a ) ( a\] p++ (t) 

+ p ++ (t) K V(\a)(a\) , (4) 

p+_(t)\a)(a\ = -ig{ae~ m + a)e m ,\a){a\}p+-{t) 

+ p + -(t)KV(\a)(a\) . (5) 

Here [•, •] ({•, •}) denote (anti-)commutator brackets. Ad- 
ditionally equations for p (t) and p |_ (t) are obtained 

by substituting g — > — g in Eqs. (4) and (5) respectively. 

These decoupled equations can be solved using the su- 
peroperator method [21, 22] assuming that the cavity is 
initiated in the coherent state \a ) 
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Next, using the commutation relations 

[a-,a f -] = 1 , [a, -a f ] = -1 , 

we obtain 

[D,A] = -A, [D,X] = X, [D,y] = 
[A,B] = 0, [X,y] = -8{l- cos St) 



-y, 
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and later use the Baker-Campbell-Hausdorff formula to 
decompose (6) into *t«] 
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Extending this treatment to two qubits in a single cav- 
ity requires taking another copy of the interaction Hamil- 
tonian (2). The only difference is that now there will 
be more Hamiltonians Vj, acting separately on different 
qubit states and jointly on the same cavity state. The 
solutions are obtained analogously to a single qubit case 
(see Appendix). 

Entanglement Evolution. Using the approach first pro- 
posed in [23], we can now quantify the degree of entangle- 
ment of a 2 x 2 system by means of concurrence defined 
as 



C = max ^0, 



0- 



where Xi are the descending eigenvalues of the real matrix 
R = {o-y ® <T y ) p* [<j y ® Gy ) p. We will assume that the 
qubit pair is initialised in either of the generalised set of 
Bell states 
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and afterward it is evolving according to the dynamics 
given by (3). As a result, the concurrence is a non-trivial 
function of time 

C = s in26le~ 8(1_cos ' 5t)(! ' 1±92)2 ( Kt+e_ ' tt ~ 1 ) / ( K2 ' 52t2 ) , (7) 

where the upper (lower) sign is used to denote the en- 
tanglement evolution of the \& (<!>) states. The graph for 
the evolution of both of these is plotted in Fig. 2. By 
approximating each of the peaks with a Gaussian we find 
that every consecutive maximum will have the form 
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where we used the standard deviation r„ to be a measure 
of every consecutive revival time given by 
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FIG. 2: Concurrence for the Bell *P (red) and Bell $ (blue) 
states. We see that the amplitude of variation is significantly 
smaller and the entanglement recovery speed are greater in 
case of the later ones. Dashed lines are Gaussians with stan- 
dard deviations given by equation (8). Plots made for g% = 1, 
<?2 = 0.5, and k — 1. 




FIG. 3: Concurrence for the Bell * (red (i) , yellow ( "' and 
pink^"') and Bell $ ( blue (i) , green c<i) and gray (i ">) states. 
Plots made for k = 1 (i) gi = 1, gi = 0.5, (ii) gi = 1, gi = 1 
and (iii) g\ = 2, g% = 0.1. If the coupling strengths are the 
same (ii) $ will experience on changes. If the relative coupling 
strength is large (iii) the concurrences for $ and 9 are very 
similar. 



Eq. (7) displays a number of striking properties. 
Firstly, after an initial sharp decrease the concurrence 
periodically recovers its initial value sin 20 never exceed- 
ing it throughout. This confirms the previous result that 
qubit-qubit entanglement enhancement is not possible 
in this system [18, 19]. Secondly, the entanglement ex- 
hibits oscillatory behaviour showing periodic revivals at 
5t = 2?i7r, with the revival time intervals r n -> oo as 
n,t — > oo. Thirdly, the greater the rate of cavity de- 
cay k (Figure 4) or the degree of detuning 5 (Figure 
5), the quicker is the recovery of the initially entangled 
state. The reason for that is that with greater k the 
cavity eigenmode field deplets quicker and so the chance 
for qubits to interact with the quantum field decreases. 
This effect is enhanced if the qubits are detuned from 



4 




15 



FIG. 4: A 3D plot of concurrence of VP as a function of time 
and detuning. For detunings S < 1 we observe the entangle- 
ment sudden death behaviour already noted by [18]. For in- 
creasing values of detuning, the concurrence function reaches 
the steady state maximum value quicker. Plots made for 
gi = 1, gi = 0.5, and k = 0.1. 

the quantum eigenmode inhibiting interaction. As a re- 
sult both of these effects lead to a decreased opportunity 
of disentanglement. Moreover, unlike Refs. [18-20], we 
have chosen to work with an arbitrary initial coherent 
state amplitude a (0) ^ to observe that its value plays 
no role in the qubit-qubit entanglement evolution, thus 
making this result universal for all cavities. 

Finally we also find, in line with Refs. [19, 20], that 
qubits initialised to different Bell states respond differ- 
ently in this system. In Refs. [19, 20] the authors claimed 
that the concurrence of the <£> type states is unaffected by 
cavity dissipation. We find that this is only true if the 
qubits are equally coupled to the cavity vacuum field. 
As a result all Bell states formed with unequally coupled 
qubits will decay and be revived depending on the values 
of S and k, however the <J> states will do it at a slower 
rate and the value of concurrence will drop to a lesser 
extent (see Fig. 2). 

To simplify calculations, we have chosen the regime 
i>w,5>3>7, These conditions can be realized 
in two types of solid-state qubits. In a superconducting 
qubit coupled to a microwave cavity [24] , one can achieve 
lo w 5 GHz, g w 100 MHz, and 7 « 1 MHz. The critical 
parameter here is the qubit dissipation 7. We do expect 
however that the condition g > 7 can be relaxed with- 
out qualitatively affecting the results. Another system is 
an NV center in diamond strongly coupled to an optical 
cavity [25], used as a spin qubit. For this realization, the 
coupling strength is the crucial parameter to observe the 
entanglement revival. 

Conclusions. We have shown that a pair of qubits de- 
tuned from the quantum field, initially show a sudden 
decrease of entanglement to later gradually recover its 
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FIG. 5: A 3D plot of concurrence of \P as a function of time 
and cavity decay rate. For decay rates s < 1 we observe 
that less entanglement is lost and its asymptotic recovery is 
quicker. Plots made for gi = 1, g2 = 0.5, and 5 — 2. 

initially entangled state. Additionally we show that no 
Bell states are strictly protected, however they can be 
better secured by skillful adjustment of the detuning pa- 
rameter. 

This study can be extended by considering the dynam- 
ics of a tripartite system composed of any Bell state and 
the coherent state of the cavity. Here one can consider 
the evolution of entanglement between any selected pair 
of subsystems and study entanglement creation between 
the cavity and the qubits during the qubit-qubit disen- 
tanglement phase. These results will be published else- 
where [26]. 

Moreover, by virtue of extendibility of this model to an 
arbitrary number of qubits as well as cavities, using this 
framework one could study multipartite entanglement 
and how, depending on the conditions and parameters 
choice, entanglement could be exchanged or transferred 
between different subsystems. Additionally, already with 
three qubits in the cavity one could try to find more dif- 
ferences in evolution between two maximally entangled 
classes: the GHZ and the W states. 
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Appendix 

Here we present the equation and solutions to Eq.(3) 
for a two qubit case. 

The density matrix of the two qubits states is labeled 
by Pij;kl,a(t), where i,j = ± and k,l = ± refer to the 
first and second qubits, respectively. Additionally in ev- 
ery entry there is a distinct coherent state density opera- 
tor \oi){a\ i -. u present. Equations for individual matrix 
components read 
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The solutions to these equations can be easily obtained 
by identifying every commutator with an A and every 
anti-commutator with a B operator. The solutions are 
again exponents of operators acting on \a )( a\ which can 



be decomposed using the Baker-Campbell-Hausdorff for- 
mula. By first using the fact that operators A and B 
commute to decompose A from |S + nDt we later follow 
the single qubit case steps to obtain 
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where the solutions to the remaining six entries of the where / and c + were defined before. Additionally, in the 
density operator are found by Hermitian conjugation and above we defined 
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These solutions can be used to extend the treatment to multipartite entanglement in this system. Tracing out 
more then one unequally coupled qubits and to study the cavity leads to a two-qubits reduced density matrix 
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where we define 
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Upon imposition of initial conditions this result can be 
used to study entanglement dynamics by means of Woot- 
ters' concurrence. 



